Reset dynamics and latching in niobium superconducting 
nanowire single-photon detectors 



Authors 

1 2 3 1 2 3 

Anthony J. Annunziata, Orlando Quaranta, ' Daniel F. Santavicca, Alessandro Casaburi, ' 

12 2 1 2 2 3 

Luigi Frunzio, ' Mikkel Ejrnaes, Michael J. Rooks, Roberto Cristiano, Sergio Pagano, ' 
Aviad Frydman, 4 Daniel E. Prober 1 

'Dept. of Applied Physics, Yale University, New Haven, CT 0651 1 
2 CNR-Istituto di Cibernetica "E. Caianiello," Pozzuoli, Italy 80078 
Dip. di Fisica "E. R. Caianiello," Universita di Salerno, Baronissi, Italy 84081 
4 Dept. of Physics, Bar Ilan University, Ramat Gan, Israel 52520 

Abstract 

We study the reset dynamics of niobium (Nb) superconducting nanowire single-photon detectors 
(SNSPDs) using experimental measurements and numerical simulations. The numerical 
simulations of the detection dynamics agree well with experimental measurements, using 
independently determined parameters in the simulations. We find that if the photon-induced 
hotspot cools too slowly, the device will latch into a dc resistive state. To avoid latching, the time 
for the hotspot to cool must be short compared to the inductive time constant that governs the 
resetting of the current in the device after hotspot formation. From simulations of the energy 
relaxation process, we find that the hotspot cooling time is determined primarily by the 
temperature-dependent electron-phonon inelastic time. Latching prevents reset and precludes 
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subsequent photon detection. Fast resetting to the superconducting state is therefore essential, 
and we demonstrate experimentally how this is achieved. 

I. Introduction 

Superconducting nanowire single-photon detectors (SNSPDs) offer high detection efficiency 
for visible and near infrared photons, with high count rates, very small timing jitter, and low dark 
count rates. 1 " 5 Typical detectors consist of a current-biased superconducting niobium (Nb), 
niobium nitride (NbN), or niobium titanium nitride (NbTiN) nanowire patterned into a meander, 
as seen in Figure 1(a). In this meander geometry, the nanowire length is proportional to the 
detection area. SNSPDs are particularly useful in applications that require high-count-rate single- 
photon detection in the near infrared, such as photon-counting communication 6 and quantum key 
distribution. 7 Development of NbN SNSPDs is most advanced. Although NbN SNSPDs offer 
higher count rates than most other near infrared single photon detectors, the count rate in large- 
area meander devices is limited by the kinetic inductance of the nanowire, which is proportional 
to the nanowire length. For a small area detector (short nanowire), the count rate can be higher, 
but very short nanowires will latch into a finite voltage state instead of self-resetting to the 
superconducting state after detecting a photon. Latching precludes practical use of a small-area 
SNSPD at high count rate. 

The count rate in a properly resetting SNSPD is set by the electrical time constant r r = L K /R L , 
where L K is the kinetic inductance of the nanowire, proportional to its length, and Rl is the load 
resistance of the readout circuit. 9 Before a photon is absorbed, the nanowire is superconducting 
with dc bias current, It, that is less than the critical current at the base temperature, I co . An 



absorbed photon will create a localized hotspot in the nanowire, which has a finite resistance, Rd, 
that drives the bias current into the load, Rl. The equivalent circuit for a device with a finite 
resistance hotspot is seen in Figure 1(b). In the desired mode of operation, the SNSPD will self- 
reset to the superconducting state. This occurs because the hotspot cools quickly after the current 
is shunted out of the nanowire, so that Rd abruptly returns to zero. The device current, Id{t), then 
exponentially returns to its initial value before photon absorption, h, with a time constant of r r . 
The equivalent circuit for this current return stage is shown in Figure 1(c). Full reset requires a 
time of approximately 3 r r so that Id(f) will be restored to approximately 95% of h. This is 
necessary in order to have a high probability of detecting the next photon that is absorbed in the 
nanowire. 10 The maximum count rate is therefore the reciprocal of this full reset time, ~ (3 T r )~ l . 
The reset time can be reduced by increasing R L 10 ' 11 or decreasing L K 12 ' 15 . If r r is too small, 
however, the device will not self-reset but instead will latch into a finite voltage state where it is 
not sensitive to photons. 10 ' n ' 16 
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FIG. 1: (a) scanning electron micrograph of a meandered Nb SNSPD with a detection area of 16 
(j,m . The dark regions are Nb, while the sapphire substrate is lighter; (b) equivalent circuit for 
the nanowire and readout with a photon-induced hotspot resistance, Rd. Prior to photon 
absorption, Rd = 0, and h flows fully through L K . The nanowire is outlined by the dotted box; the 
arrows indicate the flow of current, (c) equivalent circuit as the bias current is returning to the 
device after the hotspot has cooled back to the zero-resistance state. For this stage, the inductive 
time constant of this circuit, LyJRl, sets the time-scale for reset of the detector, which 
characterizes the decay of the load current toward zero. 

In this article, we study the reset and latching dynamics of Nb SNSPDs using experiments 
and numerical simulations. We compare these results to our measurements of NbN SNSPDs, and 
to those of other groups. 1 " 4 ' 9-11 Reset and latching have been studied in NbN SNSPDs. 11 The 
present study is the first to examine the reset and latching dynamics in Nb SNSPDs. 11 A Nb 
SNSPD has less kinetic inductance than a NbN SNSPD of the same geometry, however, Nb has 
a longer electron-phonon time than NbN. As will be shown, the kinetic inductance and the 
electron-phonon time both play a significant role in determining the reset dynamics, which 
makes this comparison instructive. This work should be relevant to understanding the limits on 
the count rates of SNSPDs made from other materials 5 ' 15 ' 16 The present work also explains our 
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previously reported experimental results for Nb SNSPDs. ' 
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II. Background 

Models for the operation of NbN SNSPDs have been presented. 10 ' 11 ' 19 " 21 Ref. 1 1 studied the 
phenomenon of latching in detail. The analysis in Ref. 1 1 uses a phenomenological model of the 
heating and cooling of the photon-created resistive hotspot to determine if a NbN SNSPD will 
latch. According to Ref. 1 1 , electrothermal feedback creates a situation where a resistive hotspot 
is either stable or unstable in the steady state, depending on I h and r r . Formation of a stable 
hotspot, known as latching, precludes operation of the SNSPD. Thus, in the case of the SNSPD, 
parameters that give "unstable" hotspot behavior are desired. 

Solutions to the model in Ref. 1 1 were obtained analytically by determining the stability of 
the hotspot under small sinusoidal perturbations. This type of small-signal stability analysis does 
not model the time-dependent formation, evolution, and decay of the hotspot. The predictions of 
the model in Ref. 1 1 were fit to data from NbN SNSPDs. By varying several of the 
phenomenological parameters of the model, good agreement between the model predictions and 
experimental data was found. However, some of the phenomenological parameters used in Ref. 
11, notably the "hotspot temperature stabilization time," are not directly connected to the 
microscopic physical processes, such as electron-phonon scattering, phonon-escape and electron 
diffusion. These physical processes govern energy relaxation in superconducting thin films, and 
are important for understanding other non-equilibrium superconducting devices such as hot 
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electron bolometers and transition edge sensors. ' Identifying which of these physical 
processes plays the key role in the SNSPD is important. Additionally, in Ref. 1 1, the values of 
some parameters that are obtained by the fitting procedure appear to be different from those of 
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independent measurements. " ' " For example, the results in Ref. 1 1 based on fitting imply a 
value of the thermal conductivity, k, in NbN of approximately 0.0017 W/K-m. Direct 
measurements of the thermal conductivity in NbN, and calculations based on the Wiedemann- 
Franz law, have obtained k ~ 0.16 W/K-m. 10 We do believe that the model in Ref. 1 1 provides 
guidance that is useful in understanding trends within the data presented, and also may capture 
physical effects occurring at the superconducting-normal interface of the hotspot that may be 
particularly important to understand latching in NbN SNSPDs, with their higher resistivity and 
shorter length normal-superconducting interfaces. 

In the present work, we develop a model to analyze latching in Nb SNSPDs based on 
microscopic physical processes. These are well-known to govern thermal relaxation in 
superconducting thin films and nanowires. We find that a photon-created hotspot can stabilize to 
either a finite resistance or cool back to zero resistance, depending on I h , L K , and R L . The factors 
that determine whether or not a hotspot will latch are identified by examining the full dynamics 
of the hotspot formation and how the dynamics depend on I h , L K , and R L . We find that almost all 
of the energy stored in the kinetic inductance of the nanowire, ViLxh , is dissipated into the 
hotspot as Joule heat when the hotspot forms. This inductive energy dissipation occurs over a 
time scale that is significantly less than 1 ns in most Nb devices. We conclude that this total 
inductive energy determines whether a device will latch for a specific value of the current return 
time, r r = L K /R L . The predictions of our model agree well with measurements of Nb SNSPDs 
with no free parameters. The parameters we use in our model are based on independent 
measurements. 

We find that Nb SNSPDs are significantly more susceptible to latching than NbN devices 
because cooling by phonon emission in Nb is much slower than cooling by phonon emission in 
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NbN. Although this slower cooling is less desirable for high count rate operation, it makes Nb a 
good material in which to study latching. The conclusions we draw from our model of Nb 
SNSPDs should be applicable to SNSPDs fabricated from other materials as well. We also tested 
NbN devices to verify that our measurement methods did not introduce any spurious behavior in 
our measurements of Nb SNSPDs. Proper shielding, filtering, and clean, low noise microwave 
design are critical to these measurements. Our measured results for NbN are similar to those 
reported in Ref. 11, confirming the validity of our measurements. 

III. Devices Tested 

The fabrication procedures for the Nb and NbN devices measured in this work are 
described in Refs. 18 and 24, respectively. All Nb nanowires were ~ 7.5 nm thick, 100 nm wide 
and approximately 110 Q/square just above the critical temperature, T c . All NbN nanowires were 
~ 5 nm thick, 130 nm wide and approximately 900 Q/square just above T c . The devices studied 
had good (~ 5%) detection efficiency for 470 nm photons and uniform line width; no significant 
constrictions were apparent. The detection efficiency of these devices can likely be improved 
significantly by using optical structures designed to increase the absorption in the Nb film, as 
was done for NbN SNSPDs in Ref. 4. A summary of the devices tested is given in Table 1. The 
variation in T c , and therefore of I co , for the Nb devices appears to be due to variation in thickness 
between samples, which affects the superconductivity strongly because the devices are very 
thin. 15 The measurement setup is described in Ref. 18 and an equivalent circuit is given in Figure 
2(a). A key feature of the experimental setup is a set of cryogenic, remote controlled RF 
switches. These enable measurement during one cooldown of a given device with various shunt 
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resistors in parallel with the transmission line but close to the device. This gives total load 
resistances Rl = 50 (with no shunt), 33, 25 (with a 50 Q shunt), 17, or 15 Q. We measured L K 
independently by incorporating each device into a resonant circuit with known capacitance and 
measuring the resonant frequency for a range of temperatures below T c and for a range of 
currents below I co using a network analyzer. The measured values at 1.7 K are reported in Table 
1 for I h = 0. The dependence of L K on I h is weak, varying by ~ 5% in Nb nanowires as h is 
increased to just below I c . This will be discussed separately in a future publication. 



Parameter 


Nb (A) 

10 x 10 |im 2 


Nb (B) 

4x4 |am 2 


Nb (C) 
10 \im line 


NbN (D) 

5x5 |im 2 


NbN (E) 
5 |imline 


length: L 


500 |im 


80 \im 


10 \xm 


105 |im 


5 |im 


T c 


4.5 K 


3.9 K 


4.0 K 


10 K 


10 K 


/ C (1.7K) 


8.2 |iA 


6.2 nA 


6.4 |iA 


26.2 |iA 


25.4 \iA 




235 nH 


60nH 


15 nH 


120 nH 


16 nH 



Table 1: Parameters of the devices studied in this work. The reported values of the kinetic 
inductance for each device in the table include approximately 10 nH of inductance from the 
measurement leads. 28 The sheet resistances are approximately 110 Q/square and 900 Q/square 
for Nb and NbN, respectively. 
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IV. Model 

The simulations we have performed are based on numerical solutions to a two-temperature 
model of heat flow similar to the model used to analyze NbN SNSPDs in Ref. 20. In two- 
dimensions, the governing equations of the model are: 



r e -ph( T e) C e (T e ) 

1 cm /„ „ x i 



L e-ph\ 1 e) K - / ph\ 1 ph) L esc 

where T e = T e (x,y,t) and T p u = T p h(x,y,t), are the electron and phonon temperatures, T Q is the 
substrate temperature, T e - P h{T e ) is the electron-phonon inelastic scattering time, r es is the escape 
time for phonons (equal to 40 ps in all simulations) , D e and D ph are the electron and phonon 
diffusivities, C e { T e ) and C p h(T ph ) are the heat capacities of the electrons and phonons per unit 
volume, jjx, y, t) is the current density, and pjx,y,t) is the resistivity. In all simulations of Nb 
SNSPDs, we use r e . ph (T e ) = r e . ph (6.5 K)(6.5 ¥JT e f with r e . ph {6.5 K) = 2.0 ns, 23 D e = 1 cm 2 /s, 27 
D ph = 0.1 cm 2 /s, C ph {T ph ) = 9.8xl0" 6 r p/) 3 (J*cm" 3 K _1 ). 10,29 For the electronic heat capacity, in the 

2 3 1 27 

normal state we use C eN { T e ) = 5.8x10" T e (J*cm" K" ), and in the superconducting state we use 
C eS (T e ) = 0.92*exp(-A(T e )/k B T e ) (J^cm^K- 1 ), 27 - 30 where A(T e ) = L76k B T*(l-T e /T c ) m . 29 

The total current flowing through the device, Id(t), is determined by the readout circuit 
(Figure 2(a)) and therefore obeys the equation: 

UR d (0 + R L )i d (0 + RS, (3). 

L K L K 
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which is obtained using Kirchhoff s laws, where Rd(t) is the resistance of the nanowire and 
where h = Vb/Rb in Figure 2(a). The spatial distribution of the current density, jjx,y,t), is 
determined by the spatially-dependent resistivity of the device, pdix,y,t). If the coordinate system 
is oriented such that the positive x-axis is along the length of the nanowire in the direction of 
current flow, the total device current at position (x) for a wire of width w, and thickness d much 
smaller than the magnetic penetration depth, is given by: 

I d (t) = I d (x,t) = d\j d (x,y,t)dy , (4) 



which, by conservation of charge, must be equal at all points (x). The local resistivity, pjx,y,t ), 
will depend on the whether the point (x,y) in the material is in the superconducting or normal 
conducting state. In our model, the resistivity is defined by: 

p d (x, y, t) = p [l - [H (T c - T e (x, y,t) ) • H (j c - j d (x, y, t))]], (5) 

where H is the Heaviside step function and p is the normal state resistivity of the film. Thus, 
Pd(x,y,f) is equal to zero or p Q , depending on temperature and current density. Since only those 
sections of the strip at point (x) that are normal for all values of (y) at (x) will contribute to Rd, it 
follows that for a wire of width w, 

aw 

where the normal length, l norm (i), is calculated numerically and is the length over which the 
resistive hotspot occupies the entire width of the wire. Typically, the maximum value of l norm is 
much less than the total nanowire length, /. Finally, we calculate the effective critical current of 
the device as a function of time, I c {t). The effective critical current is defined as the minimum 
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critical current along the length of the nanowire using the Ginzburg-Landau expression for the 
temperature dependence: 



3/2 



I c (t) = min 



w / 

, djj c (T e =0) 1 



(7) 



o V 



J 



with j c (T e = 0) determined for each device by equating the Ginzburg-Landau expression for 
I c (\.lYL)=I c (T e = Q)(\-\.1/T c f 2 to measurements of IJA.l K) and setting j c = IJwd. 29 As defined, 
the effective critical current is assumed to be independent of Id. Since the effective critical 
current is only determined by temperature, once the hotspot begins to cool below the critical 
temperature, I c {t) becomes a measure of the thermal relaxation of the highest temperature region 
of the hotspot. Thus, the time scale over which the critical current returns to near its equilibrium 
value is the hotspot thermal relaxation time, t c . This is the time required for hot electrons to 
return to near their equilibrium temperature, T . 

3 1 

We have implemented a numerical solution to these equations using MATLAB. The device 
is represented by a two dimensional grid with longitudinal grid spacing Ax and transverse grid 
spacing Ay. At each grid point, the electron and phonon temperatures are defined. From these 
temperatures, all temperature-dependent quantities are defined. When a volume-dependent 
quantity such as the heat capacity is calculated, it is calculated over the volume of the cell 
centered at the grid point (x,y) where the volume of the cell is equal to dAxAy with typical grid 
spacing of Ax = Ay = 5 nm. The absorption of a photon is simulated by increasing the 
temperature of one grid point in one time step At such that 



T e (x ,y ) = T + 



hf 



(8) 



C.(r )-(AxAyrf)' 
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where (x ,y ) is the grid point where the photon is absorbed, h is Planck's constant, and /is the 
frequency of the photon. 

V. Results 

The detection cycle in an SNSPD, as illustrated by the simulations, can be divided into 
three distinct stages, labeled on the lower curve in Figure 2(b): (i) The nanowire is biased in the 
superconducting state with a dc bias current, I h , below the critical current at the base temperature, 
ho = I C (T )', here, R d = 0. (ii) A photon creates a resistive hotspot whose resistance, Rd(t), 
increases quickly due to the fast dissipation of the inductive energy stored in L K . As a result, 
most of lb is shunted into R L , which has much lower resistance than R d ; (iii) In a self-resetting 
device, after most of the bias current has transferred to R L , the hotspot quickly returns to the 
zero-resistance state and the current slowly begins to transfer back into the device with a return 
time constant r r = RJL K . In a latching device, stages (i) and (ii) are identical to those in a self- 
resetting device, however stage (iii) does not occur and Rd remains finite. Latching is seen in 
Figure 2(b) for h = 8.1 uA. As will be shown next, latching occurs in this device (L K = 235 nH, 
R L = 50 Q) at the larger value of h because of the greater heating that occurs due to the larger 
amount of stored inductive energy, ~ ViLxh , that is dissipated into the hotspot. If the hotspot 
does not cool quickly enough as the current begins to return to the device, the device will remain 
in (latch into) the resistive state. The measured latching pulse in Figure 2(b) is for a single shot 
measurement. It has more noise than the measured self-resetting pulse because the self-resetting 
waveform displayed is an average of many pulses. The measured slow decay of h(t) in the 
latching case for t > 5 ns is due to the ac coupling of the amplifier. This slow decay is not 
observed in the simulation, which assume a dc-coupled amplifier. 

12 




FIG. 2 (color online): (a) equivalent circuit; R h » R L , R d ; (b) measured (solid lines) and 
simulated (dashed lines) output pulses, h(t), with the self-resetting case labeled by the three 
regimes of operation: (i) the device is in equilibrium with Rd = 0; (ii) initial heating: a photon has 
been absorbed, the hotspot is growing and the current is transferring into the load; (iii) the 
hotspot resistance has returned to zero, and the current is returning to the device with a time 
constant r r = L K IR L ; (c) numerical results for I c (t) (dashed lines) and Id(t) = h - /z(t) (solid lines) 
for lb = 5.0 (red) and 8.1 uA (green); the latter shows the latching case. We find hatch = 7.1 (iA 
from these simulations. We note that I c as defined here can be finite even with Rd> because I c 
is calculated according to equation (7) and therefore only depends on the electron temperature, 
and not on Id. Rd is calculated from the simulation according to equation (6). 

20 

It is essential that h be close to I co for high detection efficiency. Thus, self-reset for It, ^ 
I c is necessary for a device to have both a high count rate and high detection efficiency. In Figure 
2(c), we model the effect of increasing I h from 5 uA to 8.1 uA by computing both the effective 
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time-dependent critical current, I c (t), and the device current, Id(t), as functions of time. The fast 
dissipation of the inductive energy, ViLnh , that occurs after the photon is absorbed drives I c 
quickly to zero (Figure 2(c)) and Rd to a large value. This large value of Rd drives Id to near zero 
because I b is shunted almost entirely into the lower resistance load, Rl. This entire process occurs 
in a very short time, < 1 ns in most devices, as observed from the rise time of the current pulses 
in Figure 2(b). The cooling time in Nb, as we will show, is much longer than the sub-nanosecond 
time to dissipate the inductive energy, ViLtdb ■ 

Self-reset occurs when the trajectories of I c (t) and Id(t) are such that Id(t) < I c (t) as the 
hotspot cools. For device (A) with R L = 50 Q, Figure 2(c) shows that self-reset occurs for I b = 5 
uA but not for I b = 8.1 uA. When I b = 8.1 uA, I d (t) remains above I c (t) as the hotspot begins to 
cool, maintaining a finite hotspot resistance. The resulting dc Joule heating (occurring after the 
initial energy dissipation) enforces a stable, finite resistance hotspot state. The resistance of this 
stable state may be less than the peak resistance of the hotspot, which occurs just after the fast 
initial dissipation of the inductive energy and before any substantial cooling has taken place. 
Note in Figure 2(c) for I b = 8.1 uA that the finite value of I c for long times indicates that the 
steady-state temperature at the center of the hotspot is less than T c ; nevertheless, the hotspot 
remains resistive because Id > I c - We define a latching current, Ii atc h, as the lowest value of I b that 
results in latching. For device (A) in Figure (1) with R L = 50 Q, the model predicts that Ii atch = 
7.1 uA. 

We now show that our numerical simulations agree with experimental measurements. In 
Figure 3 we plot the normalized latching current measured for several Nb SNSPDs versus z> = 
L/c/Rl- We also plot predictions for the latching current obtained by simulating devices with the 
same parameters as those measured. In addition to the measurements and simulations for Nb 
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devices, we also plot experimental measurements of the latching current in NbN devices, to 
provide a comparison of the two different materials. Each device has a different value of L K (see 
Table 1). We vary R L to change r r in both simulations and experiments. For Nb device (A), we 

2 2 

plot predictions for two values of D e (1.0 cm /s and 0.25 cm /s) to show that faster diffusion of 
hot electrons (larger D e ) makes a device less susceptible to latching. A value of D e = 1.0 cm /s 
was used for all other Nb simulations. The simulation predictions for Nb SNSPDs and the 
experimental measurements of Nb SNSPDs are in approximate quantitative agreement. In the 
simulations, there are no free fitting parameters used. We believe the discrepancies between 
simulated and measured values of li atc h are primarily due to uncertainty in the values of the 
material parameters used as inputs to the model. ' ' 

We find in experiments that a Nb SNSPD is much more sensitive to latching than a NbN 
SNSPD with the same area and value of R L . Still, we observe that NbN SNSPDs are also 
susceptible to latching, as discussed in Ref. 1 1 . A trend we observe for both Nb and NbN devices 
is that with larger kinetic inductance, a larger value of x r is required to achieve the same 
normalized value of hatch- For larger L K , the dissipated energy ~ ViLfch is larger. For example, in 
Figure 3 the circled data points for devices (A) and (B) are both at li atc h = 0.95/ c ; device (A) has 
L K = 235 nH, R L = 25 Q, and r r = 9.4 ns while device (B) has L K = 60 nH, R L =\7 Q, and r r = 3.5 
ns. 
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FIG 3: Normalized latching current predicted by the model (open symbols) and measured (solid 
symbols) versus r r = L k IRl, obtained by varying R L in the simulations and experiments. For the 
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10x10 um meander, we plot predictions for two values of the diffusion constant: D e = 1 cm/s, 

2 2 

as used in all other simulations, and D e = 0.25 cm /s, which is used only here. 



From simulations, we find that latching occurs when the current return time, r r , is 
significantly shorter than the hotspot cooling time, r c . We define r c as the time interval between 
when I c (t) first becomes finite and when I c (t) = 0.63I CO (see Figure 3). This criterion allows us to 
compare r r to r c . Simulations show that r c depends on material parameters (e.g. the electron- 
phonon time) as well as on the inductive energy ViL K I b . After photon absorption and the fast 
dissipation of the inductive energy, the hotspot temperature is significantly greater than the 
equilibrium temperature, T . The maximum temperature is largest for the largest stored energy. 
We find from simulations of Nb SNSPDs that the temperature-dependent electron-phonon 
interaction time in Nb is the dominant component of r c , but not the only component. For small 

32 2 

temperature excursions, r c would be equal to the electron-phonon time, T e . p h(T e ) ~ (T e )~ . 



16 



However, for the large temperature changes that occur in an SNSPD, the thermal relaxation has a 
more complex behavior than simple exponential relaxation and also depends on Rj(f) and on the 
out diffusion of hot electrons. Thus, simulations are needed to predict the exact trajectories of 
I c (t) and Id(t), and therefore if a device will latch or self-reset. 

Simulations confirm that the inductive energy dissipation in the hotspot occurs on a time 
scale much shorter than r c and can be much larger than the energy of the photon that initiates the 
hotspot. From simulations, we find that the exact energy dissipated in the hotspot as it is forming 

2 2 

is equal to = l ALg(Ib -I m i„ ) + hf where I min is the minimum value of Id and hf is the energy of the 
absorbed photon, where h is Planck's constant and / is the frequency of the photon. For most 
practical devices (large area and therefore large L K ), hf « V2L K I h . In all measurements and 
simulations, we find that I min « I h , so that the total initial dissipated energy is approximately 
equal to the inductive energy, V2L K I h for all but the shortest nanowires, e.g., device (C), where 
the photon energy, hf is significant. In a self-resetting device, the total amount of energy 
dissipated in a detection event is this initially dissipated energy, ViLxh , since the hotspot quickly 
cools back to near T and Rd = after the inductive energy is dissipated. In a latching device, Id(t) 
remains above I c (t) as the hotspot begins to cool. In this case, the hotspot remains resistive as Id 
begins to increase, causing additional Joule heating after the inductive energy is dissipated. This 
additional Joule heating enforces a stable dc hotspot resistance. 

In all devices, a larger value of ViLjdb leads to a hotspot with a greater peak temperature. 
This increases the cooling time, r c , since more energy must be transferred to the phonon system 
and ultimately to the substrate in order for the hotspot to cool to near T a . Since ViLxh and r c are 
both independent of R L , we can reduce R L in order to increase z> without affecting the cooling 
time, t c . We can therefore eliminate latching in any device by sufficiently reducing R L . However, 
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reducing R L decreases the count rate, proportional to (3 r r ) _1 and reduces the output signal, ~ IbR L . 

In Figure 4, we show the effects of changing L K and lb for four different situations. We 
plot simulations of Rd(t) and I c (t) for Nb devices with R L = 25 Q. (Note that in Figure 2, R L = 50 
Q was used.) The inductive energy is varied by changing L K (compare device B with L K = 60 nH 
to device C with L K = 15 nH) and by changing I h (compare device (A) with 4 = 5 uA to device 
(A) with h = 8.19 uA). The plots of the simulated data clearly show that the peak value of Rd as 
well as the value of the cooling time, r c , increases for larger values of VJ^idb • For example, in 
Figure 4, r c increases from 2.2 to 5.0 ns when the value of V^Ljdb is increased from 5 eV for 
device (C) to 85 eV for device (A) with I h = 8.19 oA; hf= 2.6 eV in these simulations. 

In addition to circuit parameters, the thermal relaxation process also depends on material 
parameters. In particular, the cooling time, t c , decreases when the electron-phonon time is 
decreased. Also, if D e is increased, z c decreases due to faster out diffusion of hot electrons (see 
device (A) in Figure 3 where two values of D e are considered). The dependence of r c on these 
material parameters likely explains the difference in the observed latching behavior between Nb 
and NbN SNSPDs (Figure 3). NbN has a much shorter electron-phonon time than Nb at all 
relevant temperatures. This gives less latching because it more than compensates for the 
slower diffusion in NbN than in Nb and for the larger value of ViL K I b in NbN SNSPDs 
compared to Nb SNSPDs of similar geometry. 
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FIG. 4 (color online): Simulation results for Nb devices (A), (B) and (C) showing the 
dependence of Rj(t) and I c (t) on h and L K \ Rl = 25 Q; z c is labeled for device (A) with h = 8.19 
uA and is equal to 5.0 ns. For device (C) r c = 2.2 ns (not labeled). 

VI. Conclusion 

We have demonstrated that our microscopic model of the electron heating and cooling 
can predict the observed reset and latching behavior in a Nb SNSPD. Our simulations and 
measurements of Nb SNSPDs show behavior qualitatively similar to our measurements of NbN 
SNSPDs and to the results reported in Ref. 1 1 for NbN SNSPDs. In Nb SNSPDs, we find that 
the fast initial heating caused by the Joule dissipation of the energy stored in the kinetic 
inductance of the nanowire determines whether a Nb device will latch or self-reset. In Nb, 
latching occurs over a wider range of operating parameters than for NbN SNSPDs. This is due to 
the much larger electron-phonon time in Nb. Our model suggests that the maximum count rate 
for a Nb SNSPD with small L K (e.g., device (C)) is (3r c>mi „ V 1 ~ 150 MHz. To achieve this count 
rate with device (C), however, a value of Rl ~ 1 Q would be required to prevent latching if It, is 
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just below I c (see Figure 3). 

We find that the intrinsic timescale for cooling of the hot electrons is the electron-phonon 
inelastic time. This sets the lower limit to the reset time of an SNSPD. Thus, reducing the kinetic 
inductance helps only insofar as it reduces the inductive energy that is dissipated as the hotspot 
forms, which reduces the cooling time. Based on this conclusion, a NbTiN SNSPD, with a 
smaller kinetic inductance than a NbN SNSPD of the same geometry, 15 may have some 
advantages over NbN SNSPDs if the electron-phonon inelastic times are similar. However, even 
if an SNSPD were fabricated from a superconducting material with a shorter electron-phonon 
time than NbN, such as MgB2, the phonon escape time, T es , probably cannot be reduced 
significantly from its value for ultra thin NbN, where r es ~ 40 ps. Recently, new geometries that 
incorporate several nanowires in parallel have been explored. 34 ' 35 The added degrees of freedom 
for the current in these parallel topologies may allow for faster reset than standard single- 
nanowire meander-type detectors of the same detection area. 

This work was supported by NSF ECS-0622035 and an NSF Graduate Research 
Fellowship to one of us (AA). 
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